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The emission of entangled liglit from planar semiconductor microcavities is studied and the en- 
tanglement properties are analyzed and quantified. Phase-matching of the intra-cavity scattering 
dynamics for multiple pump beams or pulses, together with the coupling to external radiation, leads 
to the emission of a manifold of entangled photon pairs. A decomposition of the emitted photons 
into two parties leads to a strong entanglement of the resulting bipartite system. For the quantifi- 
cation of the entanglement, the Schmidt number of the system is determined by the construction of 
Schmidt number witnesses. It is analyzed to which extend the resources of the originally strongly 
entangled light field are diminished by dephasing in propagation channels. 
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I. INTRODUCTION 

The interaction of light and matter is a fundamental 
issue connecting elements of quantum optics and solid 
state physics. It offers a wide range of quantum effects, 
for example, the emission of various kinds of nonclassi- 
cal light. These phenomena sensitively depend on the 
interaction of light with the fundamental excitations of 
the crystal. The nonclassical correlations of such systems 
can be used for various applications, such as quantum in- 
formation processing, quantum metrology, and quantum 
communication, see e.g. [J, 

One of the most prominent quantum phenomenon 
is entanglement. It has been studied since the very 
first ideas of non-local superpositions of wave functions 
arose 0, 0] . Entanglement has been used to perform a 
number of classically impossible operations in theory and 
experiment, such as, quantum teleportation, secure com- 
munication, and distillation protocols [l|, . The latter 
ones require copies of entangled mixed states to distill 
pure entangled states 0, Q , namely Bell states . One 
problem is the feasibility of appropriate quantum mem- 
ories to store and manipulate the individual, entangled 
copies [1,0. Despite this, the determination of entan- 
glement of, in general, mixed quantum states is still a 
challenging task. 

Typically, quantum correlations are determined from 
measurements of correlation functions. Here, we aim to 
quantify the measured correlations in terms of entangle- 
ment. In the past several entanglement measures have 
been introduced jTo| . In the first instance, it is conve- 
nient to quantify entanglement for pure states only. One 
example is the Schmidt number (SN) (ll| - [l3j . For pure 
entangled states the SN counts the number of required 
superpositions of local product states to express the given 
state. A generalization to mixed quantum states can be 
achieved by a convex roof construction The SN of 
a general quantum state can be determined by making 
use of the method of SN witnesses [ll|, Eli- Recently, 



an approach based on generalized eigenvalue equations — 
so-called SN eigenvalue equations — led to a general con- 
struction scheme for SN witnesses [l6j . 

Common approaches for the generation of bipartite en- 
tangled states consider type-II parametric down conver- 
sion 17 1 or biexciton decay in quantum dots [H, [l^ • An- 



other prominent example is based on parametric phenom- 
ena in two dimensional semiconductor microcavities 
[23| . They are known to realize a strong coupling between 
cavity photons and excitons 21| resulting in an anticross- 
ing of the mixed exciton-photon modes, called lower and 
upper polariton branches. The ground state of the polari- 
tons has been studied with respect to general quantum 
properties [2^ and entanglement [13]. Stimulated scat- 
tering processes of polaritons within the lower branch 
have been shown to result in a large angle-resonant am- 
plification of the pump field [25l . [2q | . It was shown, that 
scattering processes involving both polariton branches 
can lead to the generation of photon pairs, which are 
entangled with respect to the branch index [S^, HI] . 

In the present work we show that semiconductor mi- 
crocavities can be used to generate strongly entangled 
photons and demonstrate how their entanglement can 
be identified. In our study, we apply different pump 
beams to the microcavity, which leads to the emission 
of a large number of entangled photon pairs. These pair 
correlations can be identified as a strong entanglement, 
if we decompose the emitted light into two ensembles of 
beams. The identification of these strong correlations is 
done using SN witnesses. We quantify the impact of a 
lossy channel on the strongly entangled systems by deter- 
mining the SN. This procedure is closely connected to the 
solution of the SN eigenvalue equations. As a result, we 
show that different degrees of dephasing require different 
kinds of witnesses to detect strong entanglement. 

We proceed as follows. In Sec. [TTl we briefly recapitu- 
late the physical description of the polariton formation in 
planar microcavities. The intra-cavity scattering dynam- 
ics leads to branch entangled polariton pairs, which will 
be discussed in Sec. IIIII The coupling of the polaritons 
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to radiation modes considered in Sec. IIVI yields strongly 
frequency-entangled photons. We verify their correla- 
tions by the use of entanglement and SN witnesses. In 
Sec. |Vl we study the dephasing due to the propagation 
of the frequency-entangled radiation through a linear dis- 
persive medium. Sec. IVII presents our conclusions. 



II. PLANAR MICROCAVITY MODEL 

In semiconductors the fundamental excitations are ex- 
citons. An exciton is an electron- hole pair with a radius 
Rx and a binding energy Eb = e^/(2eRx), with e being 
the static dielectric constant of the crystal. Frequently, 
excitons are approximated to be bosonic particles. How- 
ever, since they are composite particles built up from 
fermions, they have an internal structure. Moreover we 
have to take into account an effective exciton-exciton in- 
teraction prj . 

HXX ^(^Eb^Y. ^k+q&I'-q^kV • (1) 
k,k',q 

In this equation b-^^ (bj^) are bosonic annihilation (cre- 
ation) operators of excitons with wave vector k and dis- 
persion i?jif(k), and A is the sample surface. Since we 
consider planar microcavities, all wave vectors in Eq. ([1]) 
shall be in-plane. As a simplification, we work in units 
where h ^ c = 1, and we assume dispersionless excitons 
Exik) = Ex. 

To couple the excitons of the crystal to in-plane cavity 
photons with dispersion 

Ec(k) = Ec{\k\) = Ec{0)V^ + {k/kof , (2) 

where fcg = Ec{0), we have to consider the exciton- 
photon interaction. The harmonic part of this interaction 
is given by 

Hxc^^rJ2^W + ^-C-, (3) 

k 

where 2Qf{ denotes the vacuum Rabi-splitting and a^, 
(aj^) are bosonic annihilation (creation) operators of the 
cavity photons with in-plane wave vector k. It leads to 
lower (j = 1) and upper (j = 2) polariton branches 

E,{k) = l{Ec{k) + ExT ^{Ec{k)~ExY + , 

(4) 

which depend on the modulus fc = |k| only. 

Fig.[l](a) schematically shows the polariton dispersions 
i?i(k) and £'2(k) (solid lines) as well as the dispersions 
Ec{k) and Ex of the cavity photons and the excitons 
(dashed lines). Note the anticrossing of the polariton 
branches, which is due to the strong coupling of exci- 
ton and cavity photon modes. The parameter 217^ is 
oftentimes called polariton splitting, since it determines 




FIG. 1. (Color online) Sketch of the considered physical 
processes. The mset (a) shows the dispersion relations of the 
excitons, cavity photons, and polaritons. Part (b) visualizes 
the interbranch polariton pair scattering. Panel (c) depicts 
the emission and propagation of the emitted entangled light. 



the distance -E2(k) — i^i(k) when the exciton and cavity 
photon modes are resonant, Ec{k) = Ex- 

The anharmonic part of the total Hamiltonian of the 
exciton-photon interaction (saturation) is [26l . [27j 

k,k',q 

where rigat — 7/{16t:Rx) is the exciton saturation den- 
sity. Together with the exciton-exciton interaction Hxx 
it gives rise to an effective polariton-polariton interaction 

k,k',qii ,j2, 

(6) 

Here the p^j, (p]^) are bosonic annihilation (creation) 
operators of polaritons in the lower or upper branch with 
in-plane wave vector k. The effective branch-dependent 
potential V|^^k?q^"'* can be calculated through a unitary 
Hopfield transformation [2^ 

fbA _ /Milk Afi2k\ fpik\ 

{ai,J {Ahm M22wJ {P2kj ^ ' 

as 

k k' q 

^-7 = 127\/ljjk+qA/lj2k'-qA/lj3kMljjk' 

^Ps A/2 ji k+q A/l j.2k> -q Af lj3 k Afl j^k' 

+ A/2j4k'A/lj3kA/lj2k'-qA/lj^k+q^ ■ (8) 

In Eq. ([5]) we have introduced the ratio of polariton split- 
ting to binding energy, pg = 2Qji/Eii. For the matrix ele- 
ments of the Hopfield-transformation one explicitly finds 



3 




-j = 1 


5 = 




--j = 2 


5=1 




5 = -1 





0.1 

k/k„ 



0.15 



0.2 



FIG. 2. (Color online) Coefficients of the Hopfield transfor- 
mation matrix in Eq. 0, for the parameters Ec{0) ~ 1.5 eV 
and 0,R = 2 meV. The solid lines correspond to Min^ and the 
dashed ones to Mi2k- The color of the curves indicates the 
values of the normalized detuning 5. For large values of |k 
the coefficient Mfj^j^ converges to one and M^jk vanishes. 



the relations 



M22k = Miik = l/-y/l + p^ 



where 



(9) 
(10) 

(11) 



Note that in difference to the relations used in Ref. (22| 
the coefhcient Afi2k is always positive. 

In Fig. [5] we show the dependence of the squared co- 
efhcicnts M^^^l ^i2k '^^ ^^'^ modulus k of the wave 
vector k for different values of the normalized detuning 

EciO)-Ex 



6 = 



20. f 



(12) 



For large values of k the coefficient Mf^^ 1 and thus 
excitons and cavity photons do not mix. The polariton 
modes are equal to the separated exciton and cavity pho- 
ton modes. For smaller k the value of Mf-^^ depends on 
the detuning 5 and the polaritons are a combination of 
excitons and cavity photons. This mixing is due to the 
strong coupling of excitons and cavity photons. 



III. BRANCH-ENTANGLED POLARITONS 

Since we are interested in the generation of entangled 
polariton pairs, we consider a situation where the pump 
laser drives coherently the upper branch at a given wave 
vector kp. This is illustrated by Fig. [1] (b) for kp = 
O.OSfcoe^, Ec{Q) = 1.5 eV, Or = 2meV and (5 = 0. The 
pumped polaritons (solid black circles) scatter into states 
belonging to different branches, ji ^ j2 (open black cir- 
cles). Entanglement arises since both paths (indicated 
by the green lines) are indistinguishable, i. e. they are 
simultaneously phase-matched. For strong pumping we 
approximately replace the annihilation operator p2k by 



its mean field value {p2kp ) and use P|i^ 
This yields a parametric Hamiltonian 

^PP = 2 S ^^^P (^kp^,k^P,q^lkp+qP2k 



k;.ttqPk,+qPlk,-q +H.C. (13) 



that approximates the polariton-polariton interaction 
Hamiltonian (jHl) for the scattering process of Fig. [T] (b). 
Note that each pair of polariton creation operators has 
a different effective potential, such that we cannot factor 
out the effective potential as done in Ref. [l^] . Addition- 
ally, there is a mean-field shift of the branch-dependent 
energy: 



E,{k)=E,{-k) 



(14) 



with 



^^q,kp ~ o Vq^kp.O "t" ""^kp.q.O + ^q,kp,kp-q 



K'"'' k 
kp,q,q-kp 



(15) 

Assuming that the scattering wave vector q fulfills the 
phase-matching condition for the considered interbranch 
polariton pair scattering process. 



-B2(i)(kp + q) + Si(2)(kp - q) = 2^2 (kp) 



(16) 



the Hamiltonian Hp^p from Eq. ([T3| applied on the vac- 
uum state I vac) generates branch-entangled pairs of po- 
laritons in the state 

"Plkp + qP2kp-q + ^P2kp + qPlkp-q) l^ac) . (17) 



Here we introduced the parameters 



a = V, 



1222 



/3 = yj2122 



1222 Y 
1222 \2 



kp ,kp ,q7 



2122 \^ 
kp,kp,qj 

2122 \2 



kp )kp .q/ 



-1/2 

-1/2 



(18a) 
(18b) 



characterizing the properties of the material. In contrast 
to Ref. [11], the state |0) in Eq. is not a BeU state (for 
7^ /3^), which is due to the inequality of the effective 
branch-dependent potentials. As usual, + /3'^ — 1 is 
the normalization of 1 0) . 

In that the polariton energy dispersions Ei^2{k) de- 
pend on k only, the phase-matching condition is fulfilled 
if 



|kp + q|2 = |kp - q|' 



(19) 



being equivalent to q ± kp. The second part of the 
phase-matching condition in Eq. (jl6p yields 

EcilK + q|) + ExilK + q|) = ^MK) ■ (20) 

The solution of this equation gives the absolute value of 
the scattering wave vector q: 



/2i?2(kp)-i?x 



EciO) 



1 



(21) 
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FIG. 3. (Color online) Magnitude of in the (J, p3)-plane 
according to Eq. (fTSbt for £0(0) = 1.5 eV, Et = lOmeV and 
kp = 0.05^063;. The phase- matching scattering wave vector q 
follows from Eqs. ^ and (pTj) . 



Because the sign of q remains unspecified, the phase- 
matching condition is fulfilled for two equivalent inter- 
branch polariton pair scattering processes. Entangled 
polaritons in the state Eq. (jl7|) appear due to the indis- 
tinguishability of these scattering channels. 

As we have mentioned above, the value of /3 influences 
the non-local character of cf. Eq. (f^. In case 
= 1/2, we have a true Bell state, and |0) is sepa- 
rable for = or = 1. In all other cases, we have an 
entangled state as a superposition of two product states. 
Such states refer to as Bell-like states. They violate a 
Bell inequality, but not maximally Q. 

In Fig. [3] we plot the value of for the phase- 
matching scattering wave vector q following from 
Eqs. ([T^ and (PT|) as a function of the normalized de- 
tuning S and the polariton splitting to binding energy 
ratio Ps- From this figure we can deduce that the state 
described by Eq. pT)) is a true Bell state only on a spe- 
cific line in the (J,ps)-plane. For values of S and Ps apart 
from this line the state of the polariton pair is an entan- 
gled Bell-like state. Since S and Ps are determined by the 
material we are in the position to tune the entanglement 
properties of the polariton pairs. For example, we might 
consider materials, where the polariton splitting is of the 
order 257 ~ 4 meV, while the exciton binding energy ap- 
proximately is Ef, ^ 10 meV. Since the ratio of the anhar- 
monic exciton-photon interaction to the exciton-exciton 
interaction, {Anflu/ (21i?b))^, is of the order of 10~^, it is 
a fairly good approximation to omit the anharmonic part 
of the exciton-photon coupling. The particular choice 
Ps = causes a simpler effective branch-dependent po- 



tential 

^^^1^ « 12Ml,,k-hqA'/l,,k'-qAfl,3kMl,,k' , (22) 

Eb 

= a' = ^ 1/2. Obviously, 

such microcavities create polariton pairs in a Bell state 
configuration. 

Another important effect results when we apply several 
pumps with different pump wave vectors kpi , kp2 , . . . to 
the microcavity. Motivated by experiments is a pump- 
pulse train, where all kp„, n = 1,2,... are aligned in 
the same direction but have different amplitudes. Then 
we have a phase- matching condition for each kp„. Ac- 
cordingly branch-entangled polariton pairs appear for 
all phase-matching scattering wave-vectors q„ following 
from Eqs. (|19p . (j21[) by inserting the respective pump vec- 
tor kp„. Then the state of N branch-entangled polariton 
pairs takes the form 

N 

IV-) = n (""Plk,„+q„;^2k,„-q„ 
n=l 

+ /3«pL,„-hq„Plk,„-q„)|vac). (23) 

The normalization n("^n + /^n) = 1 of this state follows 
from the property + /S'^ ~ 1 for each n. 

IV. ENTANGLEMENT OF EMITTED LIGHT 

A. Frequency-entangled photons 

In the following, we consider the emission of entangled 
light from the microcavity. As shown in Ref. [l^], the 
coupling of the intra-cavity polaritons to an external field 
can be described by the quasimode Hamiltonian 

^FP = E / dwff(w)|Af,2kpaL.kP,k + H.c. (24) 

j,k 

with a frequency-dependent coupling g(a;). The creation 
operator describes an emitted photon with frequency 
Lo and in-plane wave vector k. The coupling of each 
branch to the external field is proportional to the pho- 
tonic fraction |Mj2k|^- If the |Mj2kP are of comparable 
magnitude, the branch-entanglement of the polaritons 
transfers to a frequency-entanglement of photon pairs in 
the state 

N 

\^) = n "^,+ + /^"4,- |vac), (25) 

ri=l 

where the multi-indices are defined as 

(n, ±) = (£;i(kp„ + q„), kp„ ± q„) , ^^^^ 
(n, ±) = (£^2(kpn + q«): kp" ± Q") ■ 



5 



Obviously Eii^2} (kp„ - q„) = £'1(2) (kp„ + q„) for scatter- 
ing wave vectors q„ fulfilling the phase-matching condi- 
tion (fTe]) . 

We now identify the entanglement of the multiple pho- 
ton pairs as strong entanglement by changing the point 
of view according to Fig. [1] (c). For this purpose, we de- 
compose the compound Hilbert space Ti of the emitted 
photons into two parties, Ti = %- ® ?^+, where the sub- 
spaces 'H± contain all photons emitted with an in-plane 
wave vector kp„ ± q„, respectively. In Fig. [1] (c), this 
yields two spatial subspaces for all photons emitted to 
the left-hand side or to the right-hand side. 

Let us describe this decomposition mathematically. 
We may introduce the states |0)„ for a photon with an en- 
ergy £'i(kp„-|-q„), and |1)„ for a photon energy E2{\^pn + 
q„). With these definitions we get ajj ^ol _|_|vac) = 

|0)„ (g) |1)„ and 4,- 4,+ |vac) = |1)„ ® |0)„. Thus, the 
state in Eq. (|25p reads 

N 

IV') = n ("n|0)«®|l)n+/3n|l)n® |0)„). (27) 

71=1 

The expansion of this product yields a sum of 2^ prod- 
uct states 

= \ii,...,iN) ® |1 -iTv) (28) 

with in £ {0, 1}. The sequence {in)n=i can be un- 
derstood binary representation of an integer m 
between and 2^ — 1, whereas the corresponding se- 
quence (1 — in)n=i gives the complement integer m = 
(2^ — 1) — m. As a result we obtain 

\^) = lrn\m,m) (29) 

m=0 

with coefficients 

N 

I'm = n [(1 ~ in)an + i-al^n] ■ (30) 
n=l 

The expression (1 — in)oin + inPn equals a„ for i„ = 
and /3„ for i„ = 1. The normalization condition reads 

2"-l N 

E 7l=n("" + /3^) = l. (31) 

m=0 n—1 

Note that in the form of Eq. (P^]) 1-0) is no longer a multi- 
partite product state, but a strongly entangled bipartite 
state. 



B. Identification of strongly entangled states 

To identify bipartite entanglement we use entangle- 
ment witnesses [l^H^, or, more specific, SN (Schmidt 



number) witnesses. For pure states the SN arises from 
the Schmidt decomposition of the state For exam- 
ple, if we consider the pure state jV'), cf. Eq. ()29p . the 
SN is the number of non-zero coefficients ^„i- Thus, the 
SN quantifies the entanglement based on the quantum 
superposition of the product states |m, m). SN witnesses 
can also be employed for mixed quantum states. 

The construction of SN witnesses is a challenging task. 
Recently, we have shown that one can use general Her- 
mitia n op erators to identify the amount of entangle- 
ment [l^. A (in general mixed) quantum state has a 
SN greater than r if and only if there exists a Hermitian 
operator L with 

{L)=tvpL> fr{L), (32) 

where 

fr{L) = sxrp {{i^Am^r) ■■ IV-r) SN r state} . (33) 

A SN witness can be constructed from {fr{L)I — L). Ob- 
viou sly. the case r = 1 is equivalent to an entanglement 
test [301 • A possible way to identify the value of the func- 
tion fr{L) is based on a generalized eigenvalue equation — 
the so-called SN eigenvalue equation — which takes the 
form 

L\^r) = gm + \x) (34) 

with \x) being a bi-orthogonal perturbation, cf. [l^. The 
value g is the SN eigenvalue and the vector \ipr) is the SN 
eigenvector. The largest SN eigenvalue is the value of the 
function fr{L) for the SN test in Eq. The case r = 1 
delivers the separability eigenvalue equations [30t , and we 
have shown, that they also apply to the identification of 
entanglement via negative quasiprobabilities [3l|, [s^] ■ 

Now, let us measure the entanglement with respect 
to the chosen decomposition of the Hilbert space "H. To 
determine the SN of the state, we consider the projection 
L = and obtain (L) = {ij\L\^p) = 1. For the 

function fr{L) we get 

fr{L) = max{7^^^H h7^^ : m, ^ mj for i ^ j}, (35) 

which is the sum of the r largest squared Schmidt co- 
efficients [l^. Due to the normalization of the state, 

J2m=o 7m = 1' the value of friL) is smaller than one, if 
there exist more than r values Pm 7^ 0. In conclusion, the 
considered pure state lip) has a SN of 2^, in the general 
case that all a„, /?„ 7^ for n = 1, . . . , A^. 

We conclude that the emitted light, which directly cor- 
responds to the cavity-internal quantum state, is strongly 
entangled. In order to generate such a state, the quantum 
superposition of local states \m,rn), is required at least 
r ~ 2^ times. These strongly entangled output verifies 
the internal quantum correlation between the branch- 
entangled polaritons inside the cavity structure. How- 
ever, in a more realistic scenario, we have imperfections 
causing a loss of quantum entanglement. For example, 
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the initially strongly entangled state could un- 

dergo a dephasing. In the limiting case of full dephasing, 
the state pdeph becomes 



Pdoph 



1 = X! 7ml"i,m)(TO,TO| 



(36) 



and contains no interferences of the form |m,m)(Z,Z| for 
I ^ m. In this scenario, the SN equals the minimum 
value one for the separable state pdcph- This means that 
this state is useless for any protocol based on entangle- 
ment. In the following, we will study the amount of en- 
tanglement in the intermediate region between no and 
full dephasing. 



V. DEPHASING 
A. Propagation through different linear media 

Dephasing is a critical source for the loss ofquan- 
tum correlations in the notion of entanglement |l6l . Is^ . 
One possible source of dephasing of the radiation field 
is the diflFerent propagation time in different media, cf. 
Fig. [1] (c). In the considered scenario the microcavity 
emits strongly entangled photons that shall be detected 
at a certain fixed distance. In the bipartite setting under 
study, the two parts of the entangled radiation field would 
in general propagate through different media. In case of 
pumping by a pulse train, already some small differences 
in the dispersive properties of the two media would lead 
to significant relative phase shifts and hence to an overall 
dephasing effect diminishing the entanglement between 
the output channels of the two transmission lines. 

Let us assume two media with dispersions given by 
ui±(k), where the index ± indicates the propagation in 
'H± , respectively. The Hamiltonian reads H^cph = H_ + 
where 

N 

H± = Y^ [^±(^i")4,±a„,±+w±(£^2«)al:^±a^,± (37) 

n=l 

with the energies i5i„(2n) = ^'i(2)(kpn ±qn)- Recall that 
energies for wave vectors kp„ ± q„ are identical for phase- 
matching scattering wave vectors q„. These Hamiltoni- 
ans are diagonal in the photon number basis, such that 



H±\m) = Ern,±\m) 



(38) 



with modified eigenvalues in the binary representation 
(«n)^Li of the integer m £ [0, 2^-1]: 



N 

Em,± = X! ^ i„)a;±(-Bi„) + i„uj±{E2 



n=l 



(39) 



It is obvious that the vacuum can be expressed in the 
same way using the dispersion relation a;vac(fc) = k. 



The time evolution of the initially emitted radiation is 

|V>(t)} :=:e-'(-f^-+^+)*|^A) 

m=0 



The spatial distances from the cavity to detectors in the 
left and right subspace are assumed to be equal. How- 
ever, the optical path lengths differ in both parties and 
depend on the frequency components of the propagat- 
ing fields. Effectively, the arrival times at the detectors 
differ for the different field components created by the 
microcavity system. This leads to the exponential fac- 
tor in Eq. (|40)) which takes into account the phase shift 
between the two parties of photons. 

To obtain the photon state measured by the detectors, 
we have to average over the different arrival times to 
account for the different optical path lengths. In prac- 
tice the resulting statistics depends on the details of the 
dispersive properties of both media representing the two 
propagation channels. Such a treatment must be based 
on an experimental analysis of the used channels, which 
is beyond the scope of the present paper. To demon- 
strate the basic principles, we simply suppose an equally 
distributed difference of the arrival times in the two chan- 
nels. This yields 



1 



h - t2 



dtm)){m\ 



t2~h 



(41) 



m,;=o 

X sine I Xmi - 



with 



TO, m){l, l\ 



'dt^{b- a)e-'"('^+'')/2 sine (x^^) , (42) 



where Xmi = {E^- + Ern,+ - - ^.+)' sinc(?/) = 
sv[\{y)/y. The state p(ti,i2) in Eq. (j4T|) represents the 
structure of the density operator of entangled light suffer- 
ing from dephasing. For Ai = t2 — ti — > oo we obtain the 
separable state pdcph- AH the correlations generated by 
the branch-entangled polaritons vanish in this extremal 
situation. 

Clearly there is no dephasing if the photons in both 
Hilbert spaces 'H± propagate through the same medium, 
^±{k) = u}(k), as it is perfectly realized in vacuum chan- 
nels. Under such conditions, the sum 

N 

Ern,± + = [^(-^1") + ^(-^2™)] (43) 

n=l 

is independent of to such that Xmi = for all to, Z S 
[0,2^ — 1]. The difference in the optical path lengths 
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only depends on the difference of the dispersion relations 
between left and right Hilbert space. Thus, without loss 
of generality, we can assume that we have a free-space 
propagation in H-^- and a linear medium in Ti- as antic- 
ipated in Fig. [1] (c). 



B. Detection of strong entanglement 

To quantify strong entanglement in the continuous- 
variable mixed state p(ti,t2) for finite At, we need to 
find a suitable test operator L. As we have seen in the 
case of pure states, the test operator should be closely 
related to the density operator in order to have a large 
value on the left-hand side of Eq. p2p . On the other 
hand, the value fr{L) should be as small as possible. To- 
gether this means that we should use a test operator L 
in the form 



L = 



2"~1 

E 

7n,l=0 



Tm.^i\m,m){l,l\ 



(44) 



with the positive semi-definite matrix of coefficients 

r - (r™.;)™" =0 ■ (45) 

As shown in Ref. [l^, in such a case we can obtain the 
function fr{L) just by determining the largest eigenvalue 
of all r X r principal submatrices of T. 

To construct a suitable test operator we consider the 

given State ^2) = I]m^;=o^'™.'l'^'"^)(^'^l- Let a;*-"' = 

[Xm )m=o t>e the fc-th eigenvector of the coefficient ma- 

trix {pm,i)m 1=0 fo'" ^ non-zero eigenvalue. Then we 
choose L to be the projector in the subspace spanned 
by the vectors 



,(fc) 



m, m). 



(46) 



This means Tm. 
case of a pure state we obtain 

(L) = Trp(<i,t2)i = l. 



Analogously to the 



(47) 



As long as the subspace given by all jx'-'^'') does not con- 
tain a SN r vector ji/'r), for the projection L holds: 



fr{L) = SUp(l/'r|i|V'r) < 1- 



(48) 



Hence we get a SN greater than r whenever fr [L) < 1 = 
(L). One advantage of the choice of L as a projection is 
the fixed mean value (L) representing the measurement 
outcome, whereas the other side of the inequality tests 
the SN property. 

In Figs. [5] and O we plot the SN of the state p(ti,t2) 
depending on At for different values of the normalized 
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FIG. 4. Amount of entanglement within the state p(ti,f2) 
depending on At and quantized by the SN. We apply three 
different pumps with wave vectors kp„ = 0.025nfcoei, n = 
1,2,3 to the cavity. The system parameters are EciO) = 
1.5 eV and Eb = lOmeV. The dispersion of the medium is 
chosen to be uj{k) = 0.5k. Based on our units h = c = 1, 
we choose a typical reference time scale to = 1 eV. The inset 
shows the behavior for weak dephasing. 
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FIG. 5. Amount of entanglement within the state p(ti,t2) 
for A'^ = 4 different pumps. The parameters are the same as 
in Fig. 21 aside from the additional kp„ for n = 4. 



detuning S and the ratio of polariton splitting to binding 
energy Figure S] shows the case, where the microcav- 
ity is pumped by three beams with different wave vectors 
aligned in the same direction. Hence, the maximal pos- 
sible SN of the emitted radiation is eight. Applying an 
additional pump, the maximal achievable amount of en- 
tanglement increases to 16, see Fig. [5] 

Both figures indicate that an increasing dephasing due 
to the increase of At, yields a decreasing SN. For At = 
the SN of the state /3(ti,t2) is equal to 2^, which is the 
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maximum value. The jumps of the value of the guaran- 
teed SN from r to r — 1 occur for values of At where the 
corresponding witness fails to identify a SN larger than 
r. For a fixed value of At the SN strongly depends on the 
properties, S andps, of the planar microcavity. A higher 
number of pump beams — and thus a higher initial SN — 
may significantly increase the range of At for which the 
state p(ti,t2) is still entangled (compare Figs. S] and [5]). 

VI. CONCLUSIONS 

We have discussed polariton scattering processes 
within planar semiconductor microcavities with a focus 
on the possible creation of entangled polariton pairs. In 
extension to previous works, we show that a polychro- 
matic pumping of the upper polariton branch, as moti- 
vated by experiments, leads to a simultaneous creation 
of multiple branch-entangled polariton pairs. The cou- 
pling of the intra-cavity scattering dynamics to an exter- 
nal field then transfers these kind of quantum correlations 
to frequency-entangled photon pairs. Since the entangle- 
ment properties of these photon pairs are determined by 
parameters of the device, the measurement of the photon 
correlations gives valuable information about the internal 
branch-entanglement within the microcavity. 

The simultaneous creation of photon pairs renders it 
possible to generate an arbitrary number of copies of en- 
tangled qubit states p, of the form p ® p ® ■ ■ ■ ® p. Such 
kind of states are desired to perform quantum opera- 
tions based on entanglement, such as quantum telepor- 
tation. Usually the generation of such states requires 
that a source of entangled states produces at each time 



a state p, which will be stored in a quantum memory to 
obtain the desired number of copies. Here, the number 
of pump beams or the spectral properties of a pump- 
pulse train determine the maximal number of simultane- 
ously available entangled qubits. By properly choosing 
the wave vectors of the pump field, one can optimize the 
Bell type correlations within one or more of those entan- 
gled qubits. Microcavities pumped with a single pulse of 
polychromatic light serve as generators of copies of en- 
tangled qubit states, making optical quantum memories 
superfluous. Decoherence due to the storage time in a 
quantum memory cannot occur. 

If desired, the multipartite pair correlations can be 
mapped to strong bipartite entanglement. The quantifi- 
cation of such correlations can be done via the determina- 
tion of the Schmidt number, which automatically quan- 
tifies the multipartite pair correlations and the branch- 
entanglement in the microcavity. From our results fol- 
lows that the Schmidt number of such an unperturbed 
system is maximal and it can be controlled by the prop- 
erties of the pump field. A dephasing channel diminishes 
this resource of entanglement. However, we showed that 
a high amount of entanglement can be guaranteed for a 
certain range of parameters. By using a higher number 
of pump beams or properly designed pump pulses, one 
may not only increase the initially available amount of 
entanglement, but also its resistance against dephasing. 
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